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The evolution of the ground-state nuclear shapes in neutron-rich Sr, Zr, and Mo isotopes, including both even-even and odd-A nuclei, 
is studied within a self-consistent mean-field approximation based on the D1S Gogny interaction. Neutron separation energies and 
charge radii are calculated and compared with available data. A correlation between a shape transition and a discontinuity in those 
observables is found microscopically. While in Sr and Zr isotopes the steep behavior observed in the isotopic dependence of the 
charge radii is a consequence of a sharp prolate-oblate transition, the smooth behavior found in Mo isotopes has its origin in an 
emergent region of triaxiality. 
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1. Introduction 

The study of the properties of unstable nuclei both theoreti- 
cally and experimentally is nowadays one of the most active and 
fruitful research lines in nuclear physics. Nuclear systems with 
very unusual N/Z ratios are proper candidates to get insight into 
the nuclear interaction and the impact that the associated dy- 
namics might have in low-energy observables like the ground- 
state deformation and derived quantities like moments of inertia 
or vibrational excitation energies [1]. The understanding of the 
properties of those nuclei also has important consequences in 
the understanding of other physical processes like stellar nucle- 
osynthesis of heavy elements. In particular, the neutron-rich Sr, 
Zr, and Mo isotopes with mass numbers A = 100- 110 are of 
special interest for various reasons. One of them is the role that 
these isotopes play in the nucleosynthesis of heavy nuclides in 
the astrophysical r process. Their masses and decay properties 
are an essential input to model the path, the isotopic abundances 
and the time scale of the r process in a reliable way ||2|]. In ad- 
dition, from the nuclear structure point of view, this region is 
characterized by a strong competition between various shapes, 
giving rise to shape instabilities that lead to coexisting nuclear 
shapes, as well as to sudden shape transitions in isotopic chains 
[3]. 

It has been shown that the ground states of Sr, Zr and Mo 
isotopes with N ranging from the magic number N = 50 up to 
N ~ 60 are weakly deformed, but they undergo a shape transi- 
tion from nearly spherical to well deformed prolate (or oblate) 
deformations as N = 60 is approached and crossed. Evidence 
for such an abrupt shape change includes 2 + lifetime measure- 
ments |4j,|5[] as well as quadrupole moments for rotational bands 
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|6[]. Signatures of triaxiality in Mo isotopes from spectroscopic 
studies of high-spin states have also been identified 0]. Heavier 
Sr and Zr (A ~ 110) isotopes display an axially symmetric well 
deformed shape. Above this region, it has been suggested that 
the N = 82 shell closure might be quenched far from stability. 
This quenching has been predicted by different models [8], but 
still weak experimental evidence has been found. 

Masses (nuclear binding energies) and charge radii are con- 
sidered among the most sensitive observables to explore the 
structural evolution in nuclei far from the valley of stability, 
where other types of spectroscopic experiments are prohibitive 
because of the low production rates and short lifetimes of these 
unstable isotopes. Discontinuities in the systematics of particu- 
lar mass regions suggest a change in the structure of the ground 
states among the configurations competing for the lowest en- 
ergy. Drastic changes in the mean square (ms) radius of the 
charge distributions in regions of transitional nuclei may also 
be indicators of structural changes related to the nuclear defor- 
mation. 



2. Theoretical framework 

In this work we study this phenomenology in the framework 
of the self-consistent Hartree-Fock-Bogoliubov (HFB) approx- 
imation based on the finite range and density dependent Gogny 
interaction [9] with the parametrization D1S BlOll . The HFB 
wave functions are expanded in a harmonic oscillator basis con- 
taining enough number of shells to achieve convergence of the 
results for all the nuclei studied. 

The quadrupole deformations at equilibrium are obtained 
self-consistently by a minimization procedure. Nevertheless, 
constrained calculations have also been performed to gener- 
ate potential energy curves (PEC) or potential energy surfaces 
(PES) to study the evolution with the number of nucleons of 
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the various shapes associated to the different minima displayed 
by the PECs or PESs. This has been accomplished preserving 
axial symmetry first, and allowing afterward for triaxiality in 
the most relevant isotopes where shape transitions take place. 
Hence, we have constrained the mean value of the quadrupole 
operators Q20 and Q22 (in the axial case only the former is con- 
strained), 



620 = ^<<1>hfb|2z 2 - x 2 - y 2 |<I>HFB> , 
V3 , , 

Q22 = — <4>HFBl* -rl<t>HFB> ■ 

In the triaxial figures we plot Q - y planes, where 
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and y is the angle defined as tany — Q22/Q20 [11]- With this 
definition an axially symmetric prolate mass distribution has 
7 = 0°, whereas the corresponding oblate one has y = 60°. The 
quadrupole deformation parameter B is defined in terms of the 
mass quadrupole moment 220 and ms radius (r 2 ), 
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For the description of odd-A nuclei, we have extended the 
HFB formalism using blocking techniques. The blocked HFB 
wave function of the odd-A system is in general given by 

l*HFB>o=#l*HFB>, (5) 

where |Ohfb) is a reference even-even HFB vacuum 
(AvI'&hfb) = 0) and is a quasiparticle creation operator. The 
index a stands for the quasiparticle quantum numbers char- 
acterizing the blocked state (angular momentum projection K 
and parity in the case of axial symmetry). The ground state 
of the odd nucleus is determined by finding the blocked state 
that minimizes the total energy. In the present study we use 
the equal filling approximation (EFA), a prescription widely 
used in mean-field calculations to treat the dynamics of odd nu- 
clei in a time-reversal invariant way. In this approximation the 
unpaired nucleon is treated in an equal footing with its time- 
reversed state by sitting half a nucleon in a given orbital and 
the other half in the time-reversed partner. This procedure has 
been recently justified microscopically, showing that the EFA 
can be described in terms of a mixed state density operator and 
the equations to be solved are a direct consequence of the vari- 
ational principle over the energy of such mixed state I1 1211 . It 
has also been shown that the EFA and the exact blocking pro- 
cedure are both strictly equivalent when the time-odd fields of 
the energy density functional are neglected [13]. Thus, EFA is 
sufficiently precise for most practical applications. 

In this work, we present results for one-neutron (S „) and two- 
neutron (Sin) separation energies, which can be easily calcu- 
lated from the binding energies BE, 



S„(Z,N) = -BE(Z,N) + BE{Z,N - 1), 
S 2 „(Z, AO = -BE(Z, N) + BE(Z, N-2) 



Charge radii and their differences are crucial quantities to 
study the systematics of the nuclear-shape changes as they can 
be measured with remarkable precision using laser spectro- 
scopic techniques. They are obtained theoretically by folding 
the proton distribution of the nucleus with the the finite size 
of the protons and the neutrons. The ms radius of the charge 
distribution in a nucleus can be expressed as 



ir 2 c ) = {r])z + ('j)p + (N/Z)(r 2 }„ + r 2 CM , 



(7) 



where (r 2 )z is the ms radius of the point proton distribution in 
the nucleus 



(r 2 p )z = 



j r 2 p p (f)dr 
J P P (r)dr 



(8) 



(r 2 ) p = 0.80 fm 2 and (r 2 )„ = -0.12 fm 2 are the ms radii of the 
charge distributions in a proton and a neutron, respectively. r 2 CM 
is a small correction due to the center of mass motion, which 



is evaluated according to Ref. 11411 . It is worth noticing that 
the most important correction to the point proton ms nuclear 
radius, coming from the proton charge distribution (r 2 ) p , van- 
ishes when isotopic differences are considered, since it does not 
involve any dependence on N. 

The variations of the charge radii in isotopic chains are re- 
lated to deformation effects. Other structure effects like pairing 
or spin-orbit couplings have been considered |T5[ 1611 . but in 
this work they are fixed by the Gogny interaction. For an axi- 
ally symmetric static quadrupole deformation B the increase of 
the charge radius with respect to the spherical value is given to 
first order by 



(9) 



<0 = <r%h|l + — B< 



(6) 



where usually (r 2 ) sp h is taken from the droplet model. Then, 
the measured radii have been used to estimate the changes in 
quadrupole deformation [17]. On the contrary, in this work we 
analyze the effect of the quadrupole deformation on the charge 
radii from a microscopic self-consistent approach. 

3. Results 

The PECs of the neutron-rich isotopes of Sr, Zr, and Mo from 
N - 50 to N - 68 have been calculated with the Gogny-DlS 
interaction. Our results fully agree with those of Ref. (18] that 
were obtained in the same framework. The main features of the 
shape evolution in the axial case can be summarized as follows: 
The isotopes with N = 50 - 54 show a sharp PEC around the 
spherical minimum that becomes rather shallow at Af = 56-58. 
Isotopes with N = 60 are already deformed with oblate and 
prolate minima very close in energy. In the case of Sr isotopes 
the ground state is prolate, for Zr isotopes both oblate and pro- 
late minima are found at about the same energy, while for Mo 
isotopes the ground state is oblate. Beyond N = 60 the shapes 
become stable and for heavier isotopes we obtain basically sim- 
ilar results to N = 60. It is also worth mentioning the incipient 
emergence of a spherical solution in the heavier isotopes. 
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Figure 1: Calculated S(r^) (a), S 2,1 (t>), and S„ (c) in Sr isotopes compared 
to experimental data from Ref. [27] for masses and from Ref. 1231 for radii. 
Results for prolate, oblate, and spherical minima are displayed with different 
symbols (see legend). Open circles correspond to ground-state results. 



The shape change at N ~ 60 has been predicted to a different 
extent by various theoretical models. Global calculations within 
the finite-range droplet model lfl9tl with single-particle states 
obtained from folded Yukawa predict prolate quadrupole defor- 
mations increasing smoothly from N = 50 up to N = 56, and 
then jumping suddenly to large deformations between N = 58 
and N = 60. At N = 62 the shapes stabilize until a transition 
to oblate shapes is predicted at N = 74 in Sr and Zr, and at 
N = 72 in Mo. Triaxial calculations within this approach were 
carried out in Ref. [20], showing that only Mo isotopes have a 
tendency to triaxiality in this mass region. 

The approach followed in Ref. flZlll . where PESs were stud- 
ied within a finite-range liquid drop model modified by shell 
corrections taken from deformed Woods-Saxon potentials, sug- 
gests an oblate-prolate shape coexistence in Sr and Zr isotopes 
from N = 60 - 70 with prolate ground states. Mo isotopes 
display in this case a soft behavior that develops triaxiality 
at N - 72, 74. Other calculations including rotational states 
in terms of the total Routhian surface, using non-axial Wood- 
Saxon potentials I22I1 . predict two coexisting prolate and oblate 
minima for 1Q6_116 Zr isotopes, where the prolate ground state 
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Figure 2: Same as in Fig. \l\ but for Zr isotopes. Experimental data for masses 
and radii are from Refs. 12711 and [29], respectively. 



becomes oblate beyond lll) Zr. The same calculations predict a 
y-soft triaxial minimum for 1()8 Mo. 

Relativistic mean-field calculations [23] also show increasing 
prolate deformations from N = 50 up to N = 60. In the region 
between N - 60 and N - 10 the deformations are stable for 
both Sr and Zr, whereas they change to oblate at N = 64 in Mo. 
Triaxial PESs calculated from Skyrme HF+BCS were already 
published in Ref. 02411 . where the shape evolution along Kr, Sr, 
Zr, and Mo isotopes was studied. Finally, large scale studies of 
quadrupole correlation energies and ms charge radii, based on 
both Skyrme-SLy4 [25] and Gogny-DlS 02611 effective interac- 
tions have been carried out recently. In summary, most of the 
theoretical nuclear models predict increasing deformations up 
to N — 60, where oblate and prolate shapes are developed and 
exhibit minima at close energies. Which of them becomes the 
ground state depends on subtle details of the calculation. At 
variance with the calculations mentioned above, in the present 
work we also consider odd-A nuclei in the framework of the 
EFA. 



In the next figures we can see the results obtained for 
Sir*) 50 -* = (r 2 c ) N - (r 2 c ) 50 in (a), for S 2n energies in (b), and 
for S„ energies in (c). Fig. Q] corresponds to Sr, Fig. |2] to 
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Figure 3: Same as in 
and radii are from Refs 



Fig. [T] but for Mo isotopes. Experimental 
:fs. l27ll and full , respectively. 



data for masses 



Zr, and Fig. [3] to Mo isotopes. In the three figures, results 
are shown as a function of the neutron number for all possi- 
ble values of N, even or odd, in the range of interest. Ex- 
perimental data have been taken from the mass measurements 
1 27] and from laser spectroscopy experiments carried out at 
ISOLDE/CERN and at the IGISOL facility of the University 
of Jyvaskyla (n, 28, 2911 . Focusing on the experimental data, a 
consistent picture emerges. Besides the abrupt decrease of S211 
and S„ at N — 50 corresponding to the shell closure, the evolu- 
tion of the Szn an d S„ along the isotopic chains shows a change 
in the tendency around N = 60 in Sr and Zr isotopes. This sug- 
gests a change in the ground-state shape of these isotopes. On 
the other hand, the chain of Mo isotopes shows a smoother be- 
havior. These observations are confirmed by laser spectroscopy 
experiments measuring the nuclear ms radius. In this case, the 
shape change in Sr and Zr is observed in the form of a sudden 
increase of the ms charge radii at N = 58 - 60. Again, the Mo 
isotopic chain [IV] displays a smooth behavior as the neutron 
number increases. 

In these figures we also see the calculated values using the 
oblate (down triangles), prolate (up triangles) and spherical 
(dots) shapes corresponding to the local minima in the PECs. 
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Figure 4: Q - y planes for 9 6 . 98 l° Sr isotopes with the Gogny-DlS interaction. 
The contour lines extend from the minimum up to 2 MeV higher in steps of 
0.25 MeV. 
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Figure 5: Same as in Fig. [4] but for 9 &< im < l02 Zv isotopes. 



We also plot by open circles the theoretical values that corre- 
spond to the ground states of the corresponding isotopes. In 
general, the measured S 2„ are reproduced reasonably well. The 
calculated shell gap at N = 50 is larger than observed, but this is 
a well known feature related to any mean-field calculation. The 
discrepancy can be reduced considering dynamical correlations 
beyond mean field JUS. Between N - 52 and N = 58, the 
S2n energies are underestimated by the calculations, while they 
are much better reproduced beyond N = 60. In our calcula- 
tions a change in the tendency is observed at N ~ 60, which is 
more pronounced in Sr and Zr isotopes and almost disappears 
in Mo isotopes. However, the shoulder is less apparent than 
the experimental one due to the underestimation of the data be- 
low N — 60. Notice that the open circles at N — 59 in $2n 
and S„ do not coincide with any specific deformation (oblate or 
prolate) because they correspond to the difference between the 
prolate N - 59 and the oblate N = 58 in S„ and N = 57 in 
Szn- In me case of S n energies, the amplitude of the odd-even 
staggering in S„ is well reproduced by the calculations indicat- 
ing the validity of our theoretical description of odd-A nuclei. 
In general, the agreement is fairly good below N = 50. Then, 
the calculations underestimate the measured S„ values between 
N = 50 and N = 60, being the net effect a displacement to 
slightly lower energies. The agreement improves substantially 
for heavier isotopes. 

The evolution of the nuclear charge radii can be seen in pan- 
els (a) of Figs. Q] [2] and [3] where the radius, relative to that of 
the N — 50 isotope, is plotted as a function of neutron number. 
Results for both even-even and odd-A isotopes are included in 
the plot. It is worth pointing out the almost negligible odd-even 
staggering observed both in the theoretical predictions and the 
experimental data. In Fig. [1] for Sr isotopes, we can see that 
the calculations follow nicely the measurements. There is a 
smooth increase of d(r%) up to N - 58, then a sudden jump 
occurs at N = 60 and again the increase is smooth for heav- 
ier isotopes. The encircled symbols indicate the shapes corre- 
sponding to the ground states obtained in our calculations. We 
can see that the lighter isotopes are spherical, then they change 
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Figure 6: Same as in Fig. E] but for I00 "' 10 Mo isotopes. 



into oblate very smoothly and at about N = 60 they become 
prolate. The observed jump corresponds to the transition from 
the oblate to the prolate shape since the oblate deformation is 
placed at /3 0.2 while the prolate one appears with a differ- 
ent magnitude (fl ~ 0.4). Fig. [2ja) shows similar results for 
Zr isotopes. Here, the spherical shapes account for the behav- 
ior of 8{r£) up to N = 56. From N = 56 up to N = 60 there 
is a smooth transition to oblate shapes that become the ground 
states and reproduce quite well the experiment. Above N - 60 
we obtain prolate ground states with radii in agreement with 
the observed jump. For heavier isotopes (beyond N = 66) we 
obtain again oblate shapes but there is no information in this re- 
gion. Finally, Fig. [3ja) for Mo isotopes shows that the lightest 
isotopes are spherical changing into oblate shapes and increas- 
ing the 8{r^) very smoothly. 

In general, we observe that the calculations from N = 54-60 
with spherical ground states underestimate the data for the nu- 
clear radius in the three isotopic chains. One should notice 
that these spherical solutions are not very sharp but shallow 
minima as N increases. Thus, a possible explanation for the 
discrepancy in 5(r%) could be that configuration mixing plays 
an important role in these isotopes and the actual ground state 
will have contributions not only from the spherical configura- 
tion, but also from neighbor deformed states that will increase 
slightly the charge radii. In the case of Sr isotopes, the sud- 
den change in <5(r t 2 ) occurs experimentally between N - 59 and 
N = 60, whereas theoretically it appears between N = 58 and 
N = 59. Similarly in Zr isotopes we obtain the change be- 
tween N = 60-61, while experimentally is observed between 
N = 59 - 60. We do not think this discrepancy is significant 
since it is related to the subtle competition between prolate and 
oblate shapes. We should notice that in the isotopes where the 
shape is changing we get practically degenerate energies for 
oblate and prolate deformations and thus, tiny changes in the 
details of the calculations can lead to a different ground state. 
In particular, if the ground state in N = 59 in Sr were oblate in- 
stead of prolate (they are separated by less than 0.5 MeV in the 
calculations) we will get agreement with the experiment. Sim- 
ilarly, in Zr isotopes, if N = 60 were prolate instead of oblate 
(separated by 0.2 MeV in the calculations) we would match the 
experimental jump. In Fig. |3ja) we also recognize a difficulty 
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Figure 7: Proton and neutron SPEs for 100 Zr (Z = 40, N = 60) as a function 
of the axial quadrupole moment Q20 (b). The Fermi level is plotted by a thick 
dotted line. 



in the reproduction of the data in the heavier isotopes, since the 
oblate shapes, which are the ground states, underestimate them. 
However, we can see that the prolate shapes, which are close 
in energy, agree with the data. To get a further insight into the 
reason for this discrepancy and for the fact that the observed 
jump at Af = 60 in Sr and Zr almost vanishes in Mo, we have 
performed triaxial calculations for the critical isotopes around 
Af = 60. Similar calculations can also be found in Ref. [18]. 
We can see in Figs. |H and |6] the Q-y plots for N =58, 60 
and 62 in Sr (Fig. g]) and Zr (Fig. EJ, as well as for N =58, 60, 
62, 64, 66 and 68 in Mo isotopes (Fig. [5J. In these figures we 
can see that in Sr and Zr isotopes the transition from oblate to 
prolate at Af = 60 is manifest, suddenly changing from defor- 
mations with Q ~ 2 b in the oblate sector to Q ~ 5 b in the 
prolate one. On the contrary, in Fig. |6]for Mo isotopes we see 
that the oblate shape at Af = 58 becomes gradually triaxial as Af 
increases. An island of triaxiality is apparent from Af = 60 up 
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to N = 68. We have calculated the charge radii corresponding 
to these triaxial configurations that become ground states and 
have added them in Fig. [3ja) with open squares. The new 5{r^) 
values for these isotopes are now very close to the axial-prolate 
values and agree very nicely with the experiment. One should 
notice that one important ingredient for the agreement achieved 
is that the location of the quadrupole g-value for the triaxial 
minima is much closer to the axial prolate minima than to the 
oblate ones, which are lower. 

To further illustrate the emergence and competition of de- 
formed configurations in this mass region, we show in Fig. [7] 
for the case of 100 Zr (Z =40, N =60), the proton and neu- 
tron single-particle energies (SPE) as functions of the axial 
quadrupole moment (220- Fermi levels are plotted with thick 
dotted lines. These diagrams help us to identify the regions of 
low level density, which favor the onset of deformation (Jahn- 
Teller effect), as well as to stress the important role of the inter- 
play between the proton ngg/2 and the neutron vhu/2 orbitals 
(Federman-Pittel effect) to generate deformed configurations 
1 1 ill . In the plot for the proton SPE we observe an energy gap 
below the Fermi level at Q20 = 5 b, which favors the onset of 
prolate deformation in Sr isotopes. On the other hand, above the 
Fermi energy, the low level density on the oblate sector favors 
oblate configurations in Mo isotopes. In the case of neutrons 
the energy gap below the Fermi level at ^20 = -2.5 b favors 
oblate shapes in lighter isotopes (N < 60), whereas the energy 
gap above the Fermi level at Q20 = 5 b favors prolate shapes 
in heavier isotopes (N > 60). Thus, these simple ideas offer a 
qualitative understanding of the various mechanisms leading to 
deformation in this mass region. 

4. Conclusions 

We have used self-consistent HFB calculations based on the 
interaction Gogny-DlS to study neutron separation energies 
and charge radii in neutron-rich Sr, Zr, and Mo isotopic chains. 
Our primary aim has been to search for signatures of structural 
changes, and more specifically, for shape transitions in these 
observables. We have found these correlations and specifically 
a remarkable connection between shape transitions and sudden 
changes in the behavior of the isotopic dependence of the nu- 
clear charge radii. The different sensitivities of S211 and 8(r%) 
can be understood because shape transitions take place in this 
mass region through isotopes characterized by shape coexis- 
tence, where the energies of the various shapes are almost de- 
generate. Thus, neutron separation energies are not particularly 
sensitive to these changes. On the other hand, the sensitivity to 
shape transitions is enhanced for charge radii, specially when 
the transition takes place suddenly between nuclear shapes at 
different absolute values of the deformation parameter. 

As compared to Sr and Zr isotopes, we have found that Mo 
isotopes exhibit a smoother increase in the charge radii with 
the number of neutrons, which is in good agreement with the 
experimental data. Triaxial degrees of freedom are required to 
get this agreement beyond N = 60. From triaxial calculations 
we have shown that Sr and Zr isotopes suffer a sharp transition 
from oblate to prolate shapes aiN - 60. On the other hand Mo 



isotopes display a smooth transition through a wide region of 
triaxiality. 

In summary, we have shown the ability of the isotopic differ- 
ences in nuclear charge radii to signal structural changes related 
to deformation. We have also demonstrated the capability of 
HFB Gogny-DlS calculations to reproduce those features and 
therefore to make predictions in unexplored regions. 
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